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Heavy Flavors 



Figure 1: Integrated luminosities based on projected run plans for (left) charm and (right) bottom flavor 
physics experiments [2, 3, 4, 5, 6]. At the time of this review, we are approaching the end of the first long 
shutdown (LSI) at the LHC. The Bellell integrated luminosity in the left panel is normalized to LHCb by 
the expected number of reconstructed benchmark decays, D + * —> D°n 1 and D° —> K 7Z + . 


1. Introduction 

This review outlines recent (since last year’s lattice conference) lattice QCD activity in heavy 
flavor physics and compares the precision of lattice and experiment for heavy flavor observables. 
This comparison is made using current experimental results and projections for the year 2020. 1 The 
combination of theory and experiment reveals several tensions between the Standard Model (SM) 
and nature. These tensions are reviewed in light of recent lattice results. After motivating lattice 
efforts in heavy flavor physics in the present section, section 2 focuses on processes that occur at 
tree-level in the SM, useful in precision determinations of CKM matrix elements. It is subdivided 
into leptonic (section 2.1) and semileptonic (section 2.2) decays. Section 3 discusses processes 
that occur via flavor changing neutral currents (FCNCs), i.e. rare decays (section 3.1) and neutral 
meson mixing (section 3.2). SM contributions to rare processes are small and there is a possibility 
of discernible new physics effects. Each subsection is further subdivided into D and B processes. 

A considerable collection of heavy-flavor experimental results have been amassed by the flavor 
factories (BaBar, Belle, CLEO, and BES), experiments at LEP and the LHC at CERN, a host of 
Fermilab experiments, and earlier efforts at DESY and SLAC. The Heavy Flavor Averaging Group 
(HFAG) has compiled, and averaged where appropriate, these numerous experimental results [1]. 
With plans to eclipse existing data sets, several flavor experiments are currently underway or are 
scheduled to begin in the near future. Projected integrated luminosities for this next wave of charm 
and bottom experiments are shown in Fig. 1. These plots foretell a data rich era for heavy flavor 
physics. To fully leverage these data, lattice results with comparable precision are needed. 

Having multiple experimental results allows us to extract a single CKM matrix element in 
multiple ways to search for inconsistencies. It also allows us to combine experimental data and 
theory input in a global fit to test the unitarity of the CKM matrix. These tests provide nontrivial 
consistency checks of the CKM paradigm of the SM. The hope, of course, is to find a process, or a 

1 2020 projections are based on expectations for the LHCb upgrade with 20 fb~ 1 , BESIII with 20 fb~ 1 , and Bellell 
with 50 ab~ 1 . 
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class of processes, in tension with the rest — pointing us in the direction of potential new physics. 
Despite extensive effort [1, 7, 8, 9] there are precious few tensions. The overall success of the 
CKM paradigm, combined with the lack of obvious non SM physics in high energy searches, has 
generated increased interest in rare processes. 


2. CKM Physics 

Processes that occur at tree-level in the SM allow a clean determination of CKM matrix el¬ 
ements. Ideally, we determine a CKM matrix element from multiple processes (e.g. V u /, from 
B — > XJv , B — > nlv, B —y tv , B s — > KIv, A f, —> plv, ...) and check for consistency. The study of 
these processes also allows us to verify unitarity of the CKM matrix by overdetermining the sides 
and angles of the unitarity triangle(s). These tests require the combination of results from theory 
and experiment. In order to take full advantage of a particularly well-measured quantity, lattice 
must calculate with comparable precision — and vice versa. This review compares the precision 
of lattice and experiment using decay constants and form factors. Experimental bands in plots of 
decay constants and form factors (Figs. 2, 3, 4, and 7) indicate only precision and are plotted using 
the central value of the corresponding lattice average. The search for tension between the SM and 
experiment is made at the level of extracted CKM matrix elements (Figs. 5, 6, and 8). 


2.1 Leptonic Decays 

Feptonic decays of a meson with quark content qiqj, are described by Eq. (2.1). This 
expression relates measured quantities (on the left hand side), what we calculate on the lattice 
(the decay constant f&), and the quantity we are after (the CKM matrix element |V) ; j). Higher 
order electroweak corrections in this, and all decays we will discuss, are either negligible or are 
perturbatively accounted for. 
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2.1.1 D( s ) leptonic decays 

Including the new BESIII result for 38 (D -F jttv) [10], the branching fractions of D [Sj meson 
decays are known to a precision of 5 — 5.5%. Fifetimes of the Dr s \ mesons are known to better 
than 1.5% and remaining quantities on the left hand side to better than 0.4% — branching fractions 
are the leading source of experimental uncertainty. In 2020, Bellell [11] and BESIII [12] project 
branching fraction measurements at the level of 2% and lifetimes to a precision of 1%. Translating 
these experimental precisions to target errors on the decay constants gives 3% today and 1% in 
2020. Fig. 2 plots decay constant values calculated from the lattice community together with bands 
indicating the corresponding experimental precision. 

Recent lattice works include a calculation by FNAF/MIFC [13] using the MIFC Nf = 2+1 + 1 
HISQ configurations with HISQ valence physical-mass light and charm quarks at 0.06, 0.09, 0.12, 
and 0.15 fm. They are also calculating [14] both D and B decay constants using MIFC Nf = 2+1 
asqtad configurations with asqtad light quarks and Fermilab charm/bottom quarks. This calcula¬ 
tion includes pions as light as 190 MeV and is being done at 0.045, 0.06, 0.09, 0.12, and 0.15 fm. 
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Figure 2: Recent lattice results are shown in magenta (FNAL/MILC [13], TWQCD [15], and ETM [16, 18]). 
Previous year’s results, and their averages, are shown in gray and are taken from FLAG [7]. The width of 
the solid blue band indicates current equivalent experimental precision [1, 10, 19] while the narrower dashed 
blue band shows the projected precision in the year 2020 [11, 12]. 
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TWQCD recently reported results [15] using two flavors of domain wall sea quarks with domain 
wall valence quarks at 0.06 fm and pion masses as light as 260 MeV. ETM recently calculated [16] 
D and B decay constants and /?° v; meson mixing parameters using two flavors of twisted mass 
sea quarks with an automatically @(a) improved Symanzik gauge action. The calculation used 
twisted mass valence light quarks with pions as light as 280 MeV, four lattice spacings ranging 
from 0.05 — 0.10 fm, and the ratio method to iterate from charm to bottom. They are extend¬ 
ing this work to include physical-mass light quarks [17]. Since the lattice conference they have 
also reported results [18] using Nf = 2+1 + 1 Iwasaki gauge fields with Wilson twisted mass sea 
quarks, three lattice spacings (0.06,0.08, and 0.09 fm), pions as light as 210 MeV, and using Wilson 
twisted mass light and Osterwalder-Seiler strange and charm valence quarks. The Southampton- 
Edinburgh-KEK Collaboration reported [20] on a pilot study of domain wall charm discretization 
effects in quenched calculations at 0.035, 0.05, 0.07, and 0.10 fm. Based on the results of this pilot 
study, RBC/UKQCD is using domain wall charm valence quarks on Nf = 2 + 1 domain wall con¬ 
figurations at 0.115, 0.085, and 0.065 fm with physical-mass light quarks. They plan to calculate a 
host of charm physics observables and have plans to extrapolate to bottom quark masses [21], 

2.1.2 leptonic decays 

The branching fraction for the leptonic decay B tv has been measured to 20 — 25% pre¬ 
cision [1] and B meson lifetimes are known to better than 0.5%. As with D leptonic decays, the 
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Figure 3: Recent lattice results are shown in magenta (Aoki et al. [23], RBC/UKQCD [24], ALPHA [25], 
and ETM [16]). Previous year’s results and their averages are shown in gray and are taken from FLAG [7]. 
RBC/UKQCD report separate results for chiral extrapolations to the B° (circle) and B f (plus). The solid 
blue band indicates current equivalent experimental precision [1] while the narrower dashed blue band gives 
the projected precision in the year 2020 [11]. 


branching fraction dominates the experimental uncertainty. Bellell projects [11] a significantly 
improved 3% measurement of the branching fraction by 2020. The corresponding target uncer¬ 
tainties on f[j are 10% now and 1.5% in 2020. Fig. 3 plots recent lattice calculations of //,> and 
current and projected target uncertainties for // ; . The fictitious B s decay constant is useful in ratios 
and in the parameterization of hadronic inputs, e.g. in the SM prediction of B s —> B [22]. 
Despite the favorable comparison of lattice results with experiment, the role of B decay constants 
in parameterizing hadronic uncertainties in other processes makes their improvement a continuing 
priority. 

Recent lattice works include a calculation by Aoki et al. [23] in the static limit. These results 
will anchor a heavy quark expansion, guided by results using the relativistic heavy quark action 
near the charm mass, and iterating to the anchor point via the ratio method. They simulate with 
Nf = 2 + 1 domain wall sea quarks and Iwasaki gauge fields, pion masses as light as 289 MeV, one- 
loop matching including &(a) effects, and at lattice spacings of 0.09 and 0.11 fm. Planned improve¬ 
ments include working with physical-mass light quarks and performing non-perturbative renormal¬ 
ization via RI-MOM. Their results are plotted in Fig. 3 with and without estimated G(\/mb) errors. 
RBC/UKQCD calculated fn using Nf = 2 + 1 domain wall sea quarks on Iwasaki gauge fields 
with domain wall valence light quarks and a nonperturbatively-tuned relativistic heavy quark action 
for the b quark [24]. Their simulations include pion masses down to 290 MeV and lattice spacings 
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of 0.09 and 0.11 fm. ALPHA recently reported [25] results of a calculation with &{a) improved 
Wilson light valence and Nf = 2 sea quarks, Wilson gauge action, and HQET treatment of the b 
quark with non-perturbative NLO improvement. They simulate at pion masses as light as 190 MeV 
and at lattice spacings of 0.05, 0.065, and 0.075 fm. The ETM calculation discussed above for D 
meson decay constants [16] also included results for /g (j) . 


2.2 Semileptonic Decays 

The semileptonic decay of parent meson 3 d into daughter meson 33 via the quark flavor¬ 
changing interaction < 7 ; —>■ qy is described by 


d38(3? —f Stlv) 24 7t 3 

dq 2 z^G 2 f |p^| 3 


\Vif\ 2 \f + \ 2 + ff 



( 2 . 2 ) 


Unlike leptonic decays, semileptonic decays have a kinematic dependence associated with the four- 
momentum q 2 earned away by the lepton-neutrino pan - . Contributions from the scalar form factor 
/o are suppressed by m 2 / q 2 and are negligible provided experimental data used in the combination 
are appropriately restricted in q 2 . Lattice simulations are typically performed in the parent meson 
rest frame with low daughter meson momenta, i.e. large q 2 . To determine form factors over the 
full kinematic range, mj < q 2 < (Mgs — Mg,) 2 , lattice results must be extrapolated to low q 2 . The 
extrapolation is significant, and is therefore more of an issue, for semilieptonic B is) decays. The 
z expansion [26] allows the extrapolation to be carried out in a model-independent way. Ref. [27] 
combines the z expansion with hard pion chiral perturbation theory [28] in an approach that may 
permit direct lattice calculation of form factors at low q 2 , a possibility being explored in [29]. 
Branching fraction measurements in q 2 bins combined with the calculated ^-dependence of form 
factors allows shape comparison and improves the extraction of |V]/j. 


2.2.1 semileptonic decays 

The semileptonic decays D —y nlv and D —y Klv (/ = e, ji \ decay to a tau is not energetically 
allowed) permit the determination of the CKM matrix elements \V c d\ and V rs \. Current and pro¬ 
jected experimental precisions for 38{D —t 7 lev) and 38{D —> Kev) give targets for form factor 
calculations shown in Fig. 4. 2 Relative to experiment, lattice has plenty of room for improvement 
in D semileptonic decays. Though there were no results in the past year, there are several efforts 
underway. A FNAL/MILC calculation is utilizing the MILC Nf = 2 + 1 asqtad ensembles, asqtad 
light and strange valence quarks, a FNAL charm quark, and four lattice spacings (0.045, 0.06, 0.09, 
and 0.12 fm). They plan to form ratios of form factors with decay constants to cancel charm 
quark discretization effects, then combine the ratios with the high-precision results of Ref. [13]. A 
separate effort by FNAL/MILC [31] uses the MILC Nf = 2 + 1 + 1 HISQ ensembles, HISQ va¬ 
lence quarks, three lattice spacings (0.06, 0.09, and 0.12 fm), and physical quark masses. ETM is 
calculating the form factors using Nf = 2 + 1 + 1 twisted mass gauge fields, three lattice spacings, 
and pion masses as light as 210 MeV [32]. 

2 D semileptonic determinations of \V c j\ and V cs have been largely limited to extrapolated values at q 2 = 0. Experi¬ 

mental results are typically recast as form factors using models for the q 2 dependence. The lattice efforts underway plan 
to calculate the form factors at ? 2 >0. It would behoove us to communicate to our experimental colleagues that future 
measurements would best be reported in a model independent way, i.e. as branching fractions or decay rates in q 2 bins. 
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Figure 4: Previous year’s results and their averages are shown in gray and are taken from Refs. [7, 30]. Solid 
blue bands indicate current equivalent experimental precision [1] while the narrower dashed blue bands give 
the projected precision in the year 2020 [11, 12]. 


The D s —y (j)lv decay provides an alternative to D —> Klv in the determination of V rs |. The 
first unquenched calculation of the axial and vector form factors for this decay, and their combi¬ 
nation with experiment to extract | V cs |, was performed by HPQCD in Ref. [33]. Given expected 
improvements at BESIII, theory errors dominate this determination of V cs \. Kanamori et al. are 
performing an exploratory calculation of D s —> rj 1 ' 1 / v using the QCDSF Nj = 2+1 stout link en¬ 
sembles, two SU (3)flavor symmetric sets of quark masses {M n = A7 () = 370 and 470 MeV), and 
a = 0.075 fm. They are calculating rj — r\' mixing angles and including disconnected contributions 
to D s —/ r]'lv [34]. 

A collection of leptonic, semileptonic, and other determinations of | V cc [| and V Ci are plotted 
in Fig. 5. At the current level of precision the different determinations of | | are in very good 
agreement and are all in agreement with the assumption of CKM unitarity. Given current and 
projected levels of experimental precision in D —> nlv, improvement in the corresponding lattice 
form factors for this decay is desirable. There is a slight tension between values of |F Ci | obtained 
from leptonic D s decays and those from D —> Klv and unitarity. Given this tension, it would be 
particularly interesting to improve upon the theory error associated with D s —y (j)lv (see Ref. [33] 
for discussion) to fully leverage upcoming BESIII results for this decay. 

2.2.2 semileptonic decays 

The only exclusive semileptonic decay currently used in the determination of Vj,/, is B -y nlv? 
The left-hand side of Eq. (2.2) is known from experiment to about 8 % while lattice has determined 

^Hopefully this will change with measurements of SS(B s —> Klv) and /, —> plv) at LHCb and/or Bellell. 
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Figure 5: Recent lattice results from FNAL/MILC [13], ETM [16, 18], TWQCD [15], and HPQCD [33] are 
shown in magenta. Previous year’s results are shown in gray and are taken from FLAG [7]. 


the squared form factor on the right-hand side with an error roughly twice as large. The current 
experimental error is expected to be cut nearly in half within the first 5 ah 1 of data at Bellell [11], 
i.e. within the first year or two of data taking. This is a rough characterization of the status of 
experiment and lattice as the analysis of B —> nlv is done over the full kinematic range of q 2 , as 
shown in the left panel of Fig. 6. This plot, taken from [7], shows the result of a simultaneous z 
expansion fit to lattice and experiment to extract | V u /, \. 

There are several ongoing lattice calculations of the B ^ nlv form factors. Preliminary results 
from FNAL/MILC [35] have errors matching current experimental precision and yield the most 
precise determination of \V u b\ to date. This calculation uses the MILC Nf = 2+1 asqtad config¬ 
urations with a FNAL b quark and asqtad light valence quarks, four lattice spacings (0.045, 0.06, 
0.09, and 0.12 fm), and pions as light as 177 MeV. Recent results from RBC/UKQCD [36] use 
Nf = 2 + 1 domain wall sea quarks and Iwasaki gauge fields, domain wall light valence quarks, 
and a non-perturbatively tuned relativistic heavy quark treatment of the b quark. They simulate 
at two lattice spacings (0.09 and 0.11 fm) and with pions as light as 289 MeV. HPQCD is calcu¬ 
lating form factors for this decay using the MILC Nf = 2 + 1 asqtad gauge fields, NRQCD b and 
HISQ light valence quarks, two lattice spacings (0.09 and 0.12 fm), and pion masses down to 190 
MeV [29]. This work is also exploring the feasibility of using the combination of hard pion chiral 
perturbation theory and the z expansion to include lattice simulation data at low values of q 2 . 

The right panel of Fig. 6 shows inclusive and exclusive semileptonic and leptonic determina¬ 
tions of\V ub \. Using the FNAL/MILC preliminary result, V u i, = 3.72(14) x 10 3 , a 2.4a persistent 
tension remains between the inclusive and exclusive semileptonic determinations. Though current 
experimental errors prevent \V u b\ from B tv from rivaling the precision of semileptonic deter¬ 
minations, Bellell 2020 projections for £${B —y tv) [11] combined with current lattice precision 
for fn suggest a leptonic \V u b\ with 2% precision in five years. 
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Figure 6: (left) FLAG’S [7] simultaneous fit of Belle B —► /r/v data (blue) and previous year’s lattice data 
(green) illustrates the current state of precision between lattice and experiment. Experimental results have 
errors roughly half the size of lattice errors, (right) Recent reported values for \ V u f, \ from FNAL/MILC [35] 
and RBC/UKQCD [36], and values of \V u b\ obtained from combination of SS(B — > tv, BABAR + Belle) [7] 
with recent calculations of fn by Aoki et al. [23], RBC/UKQCD [24], ALPHA [25], and ETM [16], are 
shown in magenta with previous year’s results in gray [7]. 


The semileptonic decay B s —>• Klv, which differs from B —> nlv only in its strange spectator 
quark, offers an alternative exclusive determination of \V u b\. Though not yet observed, a mea¬ 
surement is underway at LHCb and there are prospects during an Y(5S) run at Bellell. There are 
published B s -^r Klv form factor results from HPQCD [27], recent results from RBC/UKQCD [36], 
and an ongoing calculation by FNAL/MILC. Because of the similarity with B —> nlv, each of these 
calculations follows closely the same collaboration’s efforts described above for B —y nlv. In ad¬ 
dition to these works, there is an ongoing effort by ALPHA [37]. 

In addition to the obvious role it plays in tests of CKM unitarity, \V c b\ is also important in 
the determination of several phenomenologically important quantities, e.g. £${B S —y jiji) [38] and 
Ek [39]. The primary decay for the exclusive determination of \V c b\ is B -+ D lv. For B ± (the B° 
decay receives an additional correction for final state interactions) the differential decay rate is 


dF(B ± ^D°*lv) 4 n 3 

dw Cjr | Vew | 2 Ml * (Mb - M,r f Vw^l 


\V cb \ 2 xH\^(w)\ 2 



(2.3) 


where the kinematics are described in terms of w = (pb/Mb) ■ (piy /Mrr), | Qi• w| includes elec- 
troweak corrections, and Y(w)|^(w)| 2 is the conventional parameterization of the form factors. 
Most work is done in the zero recoil (D* at rest) limit where w —> 1, x( w ) L and there are 
additional simplifications (see, e.g., Ref. [40]). The left-hand side of Eq. (2.3), extrapolated to zero 
recoil, is currently know from experiment to a precision of 2.5% [1] with a 2020 Bellell projection 
of 2% for B -4 D* tv [11]. These precisions translate to target uncertainties on the lattice calcula¬ 
tion of J^(l) of 1.3% now and 1% in 2020. Fig. 7 compares lattice calculations of ^(1) to these 
target uncertainties. 
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Figure 7: Comparison of form factors at zero recoil from various lattice calculations. Recent results from 
FNAL/MILC [40, 49] are shown in magenta. Previous year’s results for B D* l v and lattice averages (gray 
bands) and are taken from FLAG [7]. Previous results for B ls \ —/ D^\lv are taken from Refs. [41, 42, 43]. 
The solid blue bands give the current equivalent experimental precision [1] while the narrower dashed blue 
bands give the projected equivalent precision in 2020 [11]. 


The decays —> D, Sj lv allow alternative determinations of | V c t,\. For B ± decay, the differ¬ 

ential decay rate is given by 


dF(B ± —y D°lv) 487T 3 

dw G 2 F \riE W \ 2 Ml(M b +M d ) 2 (w 2 - l) 3 / 2 


\v cb \ 2 \&{w )\ 2 + 


(2.4) 


Experiment has determined the left-hand side to 7.2% and Bellell projections indicate a 3% mea¬ 
surement for B — > D tv by 2020. 4 The equivalent target precisions for (1)|, 3.6% and 1.5%, are 
plotted in Fig. 7 along with existing lattice calculations. As with B —y D*lv most work to date has 
focused on zero recoil. 

FNAL/MILC published results [40] for B —y D' lv at zero recoil, obtaining an error commen¬ 
surate with the current experimental precision. Using the MILC Nf = 2 + 1 asqtad ensembles 
with asqtad light and FNAL charm and bottom valence quarks, they simulate at five lattice spac- 
ings (0.045, 0.06, 0.09, 0.12, and 0.15 fm) and with pions as light as 174 MeV, and found their 
leading error to be heavy-quark discretization effects. In an effort to reduce heavy-quark dis¬ 
cretization effects, FNAL/MILC and SWME are employing the Oktay-Kronfeld action [44] and 
have found improvements in the B meson dispersion relation and hyperfine splitting [45]. In a re¬ 
lated effort, SWME [46] plans to calculate B —> D*lv at zero recoil using the MILC Nf = 2+1 + 1 

4 Belle is expected to report an updated measurement of B —> Z)0) tv in the coming months. The B s —> D s Iv decay 
has not yet been observed. 
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Figure 8: (left) Comparison of V c f, values obtained using recent results from FNAL/MILC [40, 49], shown 
in magenta, and previous year’s results in gray [7]. (right) BABAR has measured a discrepancy of just over 
3cr in the correlated ratios of branching fractions [50] defined in Eq. (2.5). 


HISQ gauge fields, physical light quark masses, HISQ light and charm and Oktay-Kronfeld bottom 
valence quarks, and an improved heavy-light current with on-shell improvement through i 
HPQCD [47] showed preliminary results for B^ —> D {sj lv from a calculation using the MILC 
Nf = 2+1 asqtad gauge fields with HISQ light and charm and NRQCD bottom valence quarks, 
two lattice spacings (0.09 and 0.12 fm), and pions as light as 260 MeV. This work includes the 
kinematic dependence of the form factors and branching fractions for B —> Dlv and B s DJv, ra¬ 
tios of which are useful in the analysis of &{B S —> jU/i) (see, e.g., Ref. [48]). Since the conference, 
FNAL/MILC has reported results for B —yDlv form factors near zero recoil [49]. 

Fig. 8 compares values of \V c b\ obtained from lattice calculations of the exclusive semilep- 
tonic decays B —> D 1 * 1 Iv and from the inclusive decay B —> X c lv [51], revealing a persistent 3a 
tension. Also show in Fig. 8, a 3a tension between experiment and the SM has been measured by 
BABAR [50] in the ratios 


&{D^) 


&(B^D { * ] tv) 
3&{B —> T>(*)juv) ’ 


(2.5) 


In light of these tensions it is important to improve upon lattice determinations of the B —>■ D { ] lv 
form factors to fully leverage the precision of experimental results. To this end, improved deter¬ 
minations of form factor shapes and measured rates reported as functions of q 2 , or other suitable 
kinematic variable, would be particularly useful. 


3. Rare Processes 

Rare processes involve FCNCs and only occur in the SM via vacuum fluctuations. This intro¬ 
duces a loop suppression factor, often accompanied by a combination of the GIM mechanism and 
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Figure 9: From Ref. [52]: the differential branching fractions for (left) B° —» A"°p u , (mid¬ 

dle) B s — > (i)Li Li , and (right) the resulting constraint on the Wilson coefficients Cg and Ci^ p , where 
Cg = C| m +C 9 np , obtained from combination with experiment. The SM value lies just inside the 2<7 contour. 


further Cabibbo suppression. This SM suppression makes such processes potentially susceptible to 
new physics effects and provides a strong motivation for their study. 


3.1 Rare Decays 

In rare decays, as with the tree-level decays discussed above, the operator product expansion 
leads to a factorization of short distance (high energy) physics and long distance (low energy) 
physics. The model-dependent short distance physics responsible for the flavor-changing interac¬ 
tions occurs at a scale of M\y (or heavier) in the SM (models of new physics) and is characterized 
by the Wilson coefficients C-'K Long distance physics associated with the hadronization of ob¬ 
served asymptotic states is contained in hadronic matrix elements of local operators. In FCNC 
b —> s decays this is evident in the effective Hamiltonian density, 


-u&b—ts 

^ff 


4 Gp 

' V5 


W£(c,o«+c'o') : 


(3.1) 


with local operators given by, e.g., 



em h 

\6n 2 


so^P m bF»\ 0 {p 


*Yn P W) b and °1() = 


\6lZ 2 


VL(R) b (y 75 1 - 


(3.2) 

Observation, coupled with lattice calculation of the matrix elements, constrains the C- b and pro¬ 
vides evidence in support of, or against, the model responsible for the flavor-changing interactions. 

The authors of Ref. [52] calculate form factors for the decays B —y K*ll, B s —y <j>ll, and 
B s -F K*ll. The daughter vector mesons in these decays are unstable. Though associated threshold 
effects are not accounted for in current analyses, recent developments [53] suggest a path for their 
inclusion. The calculations of Ref. [52] use the MILC Nf = 2 + 1 asqtad ensembles at two lattice 
spacings (0.09 and 0.12 fm), asqtad light and strange and NRQCD bottom valence quarks, and 
pion masses down to 313 MeV. Fig. 9 reveals a slight tension between the SM and experiment, 
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and shows the resultant constraints obtained on new physics Wilson coefficients. FNAL/MILC is 
calculating form factors for B -F Kll decays using the MILC Nf = 2 + 1 asqtad ensembles at four 
lattice spacings (0.045, 0.06, 0.09, and 0.12 fm), with asqtad light and strange and FNAL bottom 
valence quarks, and pions as light as 174 MeV. As extensions of R —s- nlv calculations discussed 
above, FNAL/MILC [35] and F1PQCD [29] are also calculating the form factors for B nil. All 
FCNC rare decays are subject to potential non-factorizable effects associated with cc resonances, 
see Ref. [52] for discussion. It is not clear at present how best to deal with these effects. 


3.2 Neutral Meson Mixing 

Expanding the off-diagonal term of the neutral meson , 
weak interaction gives 


/o 


mass matrix to second order in the 


44i2 - ^ri2 — Y, 

z xjk 


,A/=1 „A/=1 
j °k 


o^ f - l \x){x\o^ f - { 


— Ex + i£ 


+L cf- 


\o7 \^°), 


(3.3) 

where Ay gives the change in flavor quantum number. The term with two Ay = 1 interactions is 
non-local for (near) onshell intermediate state \X). The calculation of these long distance effects is 
notoriously difficult, though recent progress has been made in kaon mixing [54]. The short distance 
term with a single Ay = 2 interaction is local and amenable to calculation. The space of all possible 
hadronic matrix elements of dimension 6 effective four quark operators is spanned by the basis 


vAf=2 




,0 


A f =2 


Ot f 


= Q a Lq a Q P Lq p , (ft 2 = Q a Lq P c/'Lq". 


A f =2 


(ff 

U 4 


= Q a Lq a Q p Rq P , and OJ' r 9 = Q a Lq p Q P Rq a , (3.4) 

where L/R = (1 ± 7s)/2 and color indices are explicit. Mixing matrix elements are historically 
parameterized in terms of bag parameters B ^ 0 , e.g., 

<^°| 6 »f = V°> = lf 2 ^M 2 ^ 0 B%. 


(3.5) 


3.2.1 D mixing 

In SM D mixing, the short distance contribution is doubly Cabibbo suppressed (by a factor 
of \VubV* b \ 2 ) an ^ the long distance contribution dominates. In the long distance mixing of the 
SM, the lighter down and strange quarks of the first two generations dominate and CP violation is 
negligible. As a result, strong constraints can be placed on CP violating observables associated with 
new physics, generally associated with short distance interactions described by the local effective 
four quark operators of Eq. (3.4). The lattice calculation of hadronic matrix elements of these 
operators is therefore phenomenologically important [55]. 

ETM has recently published the first unquenched calculation of the hadronic contribution to 
short-distance D mixing [56]. Using the ETM Nf = 2 configurations at four lattice spacings (0.05, 
0.07, 0.09, and 0.10 fm), ff{a) improved Ostwerwalder-Seiler valence quarks, and pions as light as 
280 MeV, they calculate the bag parameters to 3 — 5% precision. This is about five time more pre¬ 
cise than current equivalent experimental precision [57, 58] and on par with 2020 expectations [59]. 
FNAL/MILC [60] presented preliminary results for the mixing matrix elements from a calculation 
using the MILC Nf = 2 + 1 asqtad configurations at four lattice spacings (0.045, 0.06, 0.09, and 
0.12 fm), asqtad light and FNAL charm valence quarks, and with pions as light as 177 MeV. 
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3.2.2 mixing 

Because mixing is dominated, in the SM and beyond, by short distance contributions, there 
are phenomenologically relevant SM quantities to be calculated in addition to the hadronic matrix 
elements (bag parameters) discussed above. These quantities include the SM oscillation frequency 
AM( S ) and the SU(3) breaking ratio q, closely related to the ratio of CKM matrix elements \V t d/V ts \, 



AM ( ,. ) =cf i - 2 (B« ) |0f i - 2 |B;> ) ) and 


B mixing oscillation frequencies have been measured to sub-percent precisions [61]. 

In addition to the B^ decay constants discussed above, in Ref. [23] Aoki et al. also calculate 
the hadronic contribution to SM Bi s \ mixing. In Ref. [62] HPQCD gives a status report on the first 
lattice calculation of B^ s) meson mixing parameters with physical light quark masses. This calcu¬ 
lation uses three lattice spacings (0.09, 0.12, and 0.15 fm) from the MILC Nf = 2+1 + 1 HISQ 
ensembles with radiatively-improved NRQCD b and HISQ light and strange valence quarks. Early 
results show impressive precision for SM mixing parameters. FNAL/MILC updated [63] their 
nearly complete calculation of mixing parameters relevant in the SM and beyond. This calculation 
uses four lattice spacings (0.045, 0.06, 0.09, and 0.12 fm) of the MILC Nf = 2 + 1 asqtad ensem¬ 
bles, asqtad light and strange and FNAL b valence quarks, with pion masses as light as 177 MeV. 
Errors of ~ 9% are expected for the matrix element A, 10— 15% for remaining 

matrix elements, and better than 2% for £. 

4. Summary 

The past year has seen significant activity in lattice heavy flavor physics. By my count there 
are 35 recent or ongoing calculations with 20 talks or posters presented at this conference. Besides 
conveying the level of lattice activity over the past year, this review compared how we are doing 
relative to current and projected experimental precision. We are, generally, poised to fully leverage 
expected experimental results, with the possible exception of D semileptonic decays where several 
works underway should improve the situation. In addition to the determination of CKM matrix 
elements, decay constants and form factors are also used to parameterize hadronic contributions to 
many phenomenological observables, making improved calculations important independent of the 
experimental results discussed in this review. In light of recent lattice results, we reviewed several 
tensions between nature and the SM, including: leptonic and semileptonic determinations of |V Ci j|, 
inclusive versus exclusive determinations of \V u b\ and \V c b\, the correlated ratios 38(D) and M(D'), 
and branching fractions for some rare B^ decays. Continued coordination with phenomenolo- 
gists and experimentalists is crucial to extracting as much as possible from each measurement and 
calculation, and to identifying useful new observables in our search for cracks in the SM. 
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